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ABSTRACT 
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The method of discrete ordinates is extended to describe the 
steady state intensity distribution of thermal radiation and the 
corresponding depth-dependent thermal structure of a plane- 
parallel semi-infinite particulate medium in radiative equilibrium 
in the presence of an outside point source of visible radiation. 
The greenhouse factor [ BCr) / F(0) ] is found to approach a 
finite limit as T— Numerical calculations show that the 
diffuse field of visible radiation generated by multiple scatter- 
ing processes is extremely effective in heating the lower regions 
of a particulate atmosphere. For slightly more realistic single 
scattering parameters than those considered in this study, it 
would appear that the large differences in temperature between 
the surface and the effective level of infrared radiation of Venus, 
as inferred from the infrared and microwave observations, can be 
mostly accounted fori Neither high surface pressures nor internal 
heat sources are required to make the heating mechanism work 
effectively. 
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I. INTRODUCTION 


It has been demonstrated (King, 1963) that under certain 
conditions the temperature of a planetary atmosphere in radiative 
equilibrium increases with optical depth, even though the sole 
source of heating decreases exponentially with optical depth. The 
major requirement for a strong greenhouse effect appears to be 
that the ratio of the infrared to visible absorption coefficients 
be fairly large at most infrared wavelengths. Frequent and fairly 
wide infrared windows are the rule for atmospheres containing 
realistic amounts of optically active gases. In general such 
windows do not exist, however, for atmospheres composed of parti- 
cles of a few microns in diameter. Exceptions to this rule are 
likely to occur in rather narrow spectral intervals just to the 
shorter wavelength side of regions of anomolous dispersion, where 
jt he real part of the refractive index of the relevant substance 
is close to unity. For this and several other reasons to be dis- 
cussed it would seem germane to investigate the role that parti- 
culate matter plays in giving rise to strong greenhouse effects. 

In order to do this in any proper sense it is mandatory to include 
anisotropic scattering effects in solving the appropriate equation 
of radiative transfer. 
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II. SOLAR HEATING 


Consider- first the diffuse field of visible radiation 
arising from sunlight multiply scattered in a particulate atmos- 
phere. Approximating the sun as a point source, the equation 
of transfer becomes (Chandrasekhar, 1960; henceforth referred to 
as R. T.) 

<t>) 






- 4) -yrj 4 

o -i 




(i) 


where the symbols have the same meaning as in R.T.. In particular, 
ttF q is the flux of direct visible radiation from the sun which 
crosses a unit area normal to the beam. 

It will be assumed throughout this study that the atmosphere 
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in question is grey with respect to both visible and infrared 
radiation, although the two wavelength regions will be considered 
as having two different greyness parameters. By this is meant 
that the scattering and absorption properties of the atmosphere are 
the same throughout each of the two spectral regions, but that the 
atmospheric properties become different in going from one to the 

other. Thus (1) is independent of wavelength. 

* __ 

As is customary, w© assume the intensity in (1) to be of 
the form 
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and the phase function for single scattering representable by 
the finite series expansion of Legendre polynomials: 
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where ® is the angle through which radiation is singly 
scattered. As is demonstrated in R.T., equation (l) separates 
into (N + 1) linearly independent equations. In particular, 
the azimuth-independent term in (2) obeys the equation 
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The solution to (4) by the method of discrete ordinates 


is Jz T a/ 
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where the cosine of the zenith angle |j is evaluated at the 2n 
discrete intervals (i = 1, n). As in R.T., the ^ymbols 

1 f 0 and are defined by the expressions: 


where 


and 
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§fc) = Z ^ (8) 
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where 



^ - -^/j / . £%) j (9) 


and where the a^’s are the usual Gaussian weights determined in 
accordance with the relation 
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An alternative expression for is given by the recursion 
relation 

r ■ •xli+o > ‘ ■ ai> 

In both (8) and (11) 

ijso = i 

t/x) = o U<o) 

The n constants of integration in (5) may be evaluated by 

imposing the n boundary conditions 

I (&;-/*/) = O U*!) ‘ . (13) 

We shall be interested solely in the atmospheric heating 
caused by the incident solar radiation. For this purpose an 
expression is required for the flux divergence of the total 
visible radiation field. In this context the azimuth-depebdent 
terms in (2) contribute nothing, and we shall neglect them. 

The expression for the net flux 'n* ^(.T) of the diffuse 



radiation field, with the aid of (9) and the definition 


becomes 




( 15 ) 



By virtue of (8) and the relation* (R.T., page 154) 

A/ 



equation (15) becomes 




or, from (11) , (12) , and (14): 


F t ® = ffa ^ e^ r +/i, % 


(17) 


C 18 ) 


where (Z)% , the albedo for single scattering in the visible, is 
superscripted with a "V" to differentiate it from the albedo 
for single scattering in the infrared. 

^‘Phe term in (16) was omitted in a previous report 

(Samuelson, 1965) , and invalidates the subsequent derivation 
for F r (r) . Equations CIO, through C20 and also equation 96 
of the report, which depend upon the derivation, are conse- 
quently incorrect. 
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The net flux IT F 0 (r) of the attenuated direct radiation 
from the sun in the direction yu 6 and at a level X is clearly 
given by 

ECt) — -M E A 'tyo 

0 r ° e (19) 

The net flux 7 rFfr) of .the total visible radiation field must 

be the sum of and 7T/~(r) . Hence, from (18) and (19), 

KV = - o-sf ) ?{i +^xe W/ '} , (20) 

where t v is the optical depth in the visible. Differentiating 
(20) with respect to we obtain for the flux divergence (apart 
from the factor IT): 


dF(r v ) 

otr v 






( 21 ) 
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III. RADIATIVE EQUILIBRIUM : 

THE TEMPERATURE PROFILE 

i We now assume (King, 1963) that the flux divergence in the 
ivisible must equal the flux divergence in the infrared; i.e., in 
:a steady state there are no sources or sinks of radiation 
internal to the atmosphere. This is what shall be meant by an 
atmosphere in radiative equilibrium. For a quantitative com- 
parison we shall need an expression analogous to (21) for the 
infrared. 

The appropriate equation of transfer for a partially thermally 
emitting, partially anisotropically scattering semi-infinite 
atmosphere is (Samuelson, 1965; henceforth referred to as R-215) 

+1 

= X(SjA - if J/-0- otf) B(r) } (22) 

where B(D is the Planck function in units of specific intensity 
and integrated over all (infrared) wavelengths, and 

~ X H * (23) 

Jl-b 

All parameters in (22) refer now to the infrared, and are not to 
be confused with previous parameters referring only to the visible. 
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In particular €$ is the albedo for single scattering in the 
infrared. 

We may obtain an expression for fi(D by multiplying through 
both sides of (22) by d|i and integrating over the interval 
(-1 s )] s + 1) , With the aid of (23) , the identity 


H 

f %(*)%&) = 


-/ 
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(24) 


and the definition for the flux integral, 

F(f) r X f/* Ify/*) > 

-/ 

a straightforward evaluation yields 



dm 

dr 
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(25) 
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and (22) becomes 
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where 


This is to be compared with equation (3) in King’s (1963) paper. 
We now require in accordance with King that 


dF(zd _ _ dFCn) 
d t dz * 


(29) 


where 2 is the vertical distance measured from some reference 
height (e.g. Z 0 ) to the height of interest. 

By definition 

v 


dr Y = -j/ 0 (i)x s di 


and 


S 

If 


(30) 


dr XR = -/i 46) J* ; 


(31) 


where is the number of particles per unit volume at a level 

v x/? 

£ and and are the effective extinction (absorption 

plus scattering) cross-sections per particle in the visible and 
infrared respectively. We shall assume that the particles are 

s/ 

spherical and homogeneous, and that the definitions of X £ and 
x/? 

Xg in terms of the Mie parameters as- given in R-215, 
section 3, are valid. 

i 

Let 

/3 = 

-v 

Then, from (30) and (31), 

- fi > 




(32) 


(33) 
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or, if 3 is assumed to be independent of z: 


'Ey - ft c zR * 


( 34 ) 


This latter assumption is of doubtful validity for realistic 
atmospheres. 

By virtue of (21) and (29) - (34) , equation (27) becomes 


/* = m t M) -ifqfy) fp t m (y)4» 




~^ C + % y 


(35) 


where T* t TK . Since = 1 the homogeneous part of (35) is 
formally equivalent to conservative scattering. Denoting the 
solution to the homogeneous part by the subscript G, we have, by 
the method of discrete ordinates in the n^ h approximation* 

(R.T., page 154): 



COhstth-tK* / 


M, e* r 



(36) 




U- yi) j 


*There is no necessity for solving both (4) and (35) in the same 
degree of approximation. It is done here merely for the sake 
of simplicity. 
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where K p ( /O = l f , ,,, n-l) are the n-T positive nonzero roots 
of the characteristic equation [cf. (8)-(J2.)j 


/ 


Z q t&tfr j) 


(Js . < 37 > 


That K=0 is also a root satisfying (37) is a characteristic common 
to all problems formally equivalent to conservative scattering. 

It is clear that the (n+1) constants ( f> =1, . .., n-l), 

M 0 , and M n cannot all be determined from conditions (13) . We shall 
find presently in obtaining the complete solution to (35) that an 
appeal to the flux integral resolves the difficulty. 

Consider the equation [cf. eq. (35)] 



where C p and ^ are constants. Assume a solution for a particular 
integral of (38) of the form 


J/yO = TjCrffie** ) <39> 

where ^v)is a function of /* yet to be determined. In the scheme 
of the n^h approximation, upon replacing by (i « +1, ..., +n) 


12 


and substituting for X ( T ) / A i) i- n (38) according to (39), we 
obtain: 


V —i 

U+Vdtffti - J Z 0 * < 40 > 


Therefore is of the form 

= 7— ; 


(41) 


where the constants are S iven by 

£<6) = -5- T a j£(“/)$lfy) (J-Oy.^/y) , W2) 

/ 

Substituting in (40) according to (41) we obtain 


(43) 


Since (43) is valid for all we obtain, after some reduction, 




(44) 
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where 


= x y_g. £^Btk , 

if V ItMp (45) 

Equation (44) yields an NxN system of equations which may then 
be solved for the N values of £j for the argument 

We may also obtain from (45) a recursion relation satisfying 
the D’s: 

^ - #'■>£#> +* %$} ; (46> 

where ^ ^ is the Kronecker £- function. Along lines analogous 
to R.T., pages 154-155, we may use (44) and (46) to obtain the 
recursion relation satisfied by the 

^ ■ ~ yU+p L ^ ' (47> 

In particular ( =■ / ) : 

Lfy) - y ; (48) 

a relation we will find useful in resolving the value of M 0 in (36) . 

It is clear that (47) provides an alternative method of 
evaluating all (X = 0, N) once the value of is known. 
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We try a relation of the form 


i&) = £&) tfi) f i(i) , (49) 

where £(j?) is unspecified for the moment. Substituting (49) into 
(44) we obtain 





hence , 


upon setting ^*0 




(51) 


Again, substituting (49) into (47) we obtain, with the aid of (11), 
the relation 


= - ZitLzZLfA)^ JL 


i+) 


w + 




/) ^ 


(52) 


where, according to (12) and (49): 

£ 6 ) = o, < 53 > 

making (52) determinate for all £ , It appears possible to express 
at least part of (51) in terms of Chandrasekhar’s H-f unctions; 
however, in the interest of conserving space, we shall not pursue 
the matter here. 
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By virtue of (41) the particular integral (39) becomes 
Tpfy/it) - IjCpQ (a e £/) 7l ). (54) 


By analogy with (38) and (54), the solution to (35) is [cf. also 


(36)] 


i(x ? ^ = 

]/% +K-*. 




r r <- 
+ /+/V ^. L fa * 



(/ :*tj.:,±n) , 


(55) 


The value of M 0 in (36) may be resolved in the following 
way. Substituting (55) into the flux integral 

F&) = 2 ) (a ~ ^ ! y " 

A 

we obtain with the aid of (45) successively the relations 
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or, by virtue of (8) and (44) : 


\ 
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( 58 ) 


Rr)=U-a + jd a 

+ ^ + z e ' er<6 ° z($)} • 

Again, by (11), (12), and (48) , since U3 0 = 1, 

Fir) = 0-z$)F.fa/1o F Z jgt^F^r. ***}) 


(59) 


or, from (20) and (34) , 


m+FC r v ) = ^(J-eOfi^/1, . 


(60) 


In a steady" state 

%(*) FF/o) = o ; (61) 

hence M 0 is zero for O , and quite generally (60) becomes 


Flr x ,) F FCr v ) =0 


C/S/o) : (62) 
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Re-writing (55) and (59) we have 


XCrrf = +*H 

r&Vu. 


and 


+ 2L£ 

I + /“;^ 




/W=0-«JE{|^« A "v/<.r.e-^} , 

from which it follows that [cf. also (21) and (29)-(34) ] 


(64) 


dm 

dr 


= -0-a:)tf{iue k ' 3T + r.t 4 #'} 


(65) 


In a fashion completely analogous to the derivation of (57) 
et seq. we may obtain the relations 

H 

* Za.ICr^j) (; = ±! } ... } ±H) 

- 1 A 



s 


( 66 ) 


Hence, from (65) and (66), expression (26) for the integrated 
Planck function becomes 

BW = (1 ft* 1 * f My, 



The desired distribution of temperature with optical depth is 
given through the relation 


B(r) = £■ rh ; 


( 68 ) 


where <r = 5.672 x 10” 12 watts/cm 2 /deg 4 . 

It should be noted in (67) that B(X) approaches a finite 
limit as r-*co. Hence, analogous to equation (25) of King's 
(1963) paper, we may define a greenhouse factor from (64) and 
(67) by the relation 


Q(co)/F(o ) 



(69) 


The remainder of this study will be devoted to the numerical 
evaluation of some of the more relevant equations derived above 

t 

for certain physically interesting situations. The results will 
be discussed in the context of the atmosphere of Venus* 
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IV. SINGLE SCATTERING 

The atmosphere of Venus appears to support a cover of clouds 
sufficiently thick to obscure completely any trace of surface 

features. In addition this atmosphere may possibly support a 

I 

large greenhouse effect (Sagan, 1962). It is of considerable 
practical importance, then, to investigate the role that clouds 
may play in the production of a greenhouse effect in the atmos- 
phere of Venus. For the purpose of a preliminary reconnaissance 
of this general problem it would seem that a knowledge of the 
vertical temperature profiles and laws of darkening for a limited 
number of physically realistic model cloudy atmospheres, computed 
in accordance with the theory outlined above, would be of con- 
siderable interest. Before proceeding with such computations, 
however, it is necessary to adopt single scattering parameters 
appropriate to the atmosphere of Venus. 

According to Sobolev (1963) the albedo for single scattering 
in the visible is very high ( ©^*.989). Although the high 
accuracy of the observational data from which Sobolev derived his 
value of ©/ has been only imperfectly demonstrated, it is felt 
that the superior quality of his. work lends considerable weight to 
the essential correctness of the magnitude of this value. For the 
purpose of constructing a preliminary set of models the value 
O? - .99 has been adopted in this study. 


■ S, 
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In addition to the value of €&£ , it is, necessary to ass um e 
values of several other parameters relevant for single spattering 
in both the visible and the infrared. To this end we shall assume 
that the particles of interest are spherical and homogeneous. 

For one such particle the size parameter ©< is given by 




(70) 


where r is the particle radius and X is the wavelength of 
radiation of interest. We shall adopt "eff ective" -wavelengths of 
Xy = 5625 X and X^- in the visible and infrared respectively, 
yielding an effective size parameter ratio for a single particle of 


oP/o?* = j (e> . 


(71) 


For the purpose of computational ease and uniqueness of 
solution it was decided to consider only square particle size 
distributions, i.e. distributions of particle sizes that are 
flat over the ranges of interest. Under the assumption of homo- 
geneous spherical particles, the Mie theory was used to compute, 
for a variety of particle size distributions and complex indices of 
refraction, the following infrared single scattering parameters: 

(1) , the albedo for single scattering, (2) , the 

effective extinction cross-section per particle, and (3) -flo(ccQ> <B) , 

the phase function for sin'gle scattering. Eighteen values of the 
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complex index of refraction were used, covering the range 
(1.1 - O.Oli £ Yl £ 1.6 - 0.50i). In practice the particle 
size distributions considered were replaced,by size parameter 
distributions over the range (0 £ cK < o( c ) . Sixteen values of 
o( c were considered for each value of 91 . 

The computed results indicated that is very 

dependent upon cx 0 , but is almost independent of H , Hence, 
one mean phase function for each value of oq> was considered 
sufficient to represent the angular distribution of singly 
scattered infrared radiation for a wide variety of possible sub- 
stances. while three values of o<« were chosen as being rep re- 

• 'I 

sentative of a reasonable range of physically interesting c.on~ 
ditions. These latter values, along with the corresponding values 
of /9 [cf. eq. (32)], , GJ* , and e*'/ are included in Table 1. 

The visible radiation single scattering parameters required to 
complete the table were calculated for an assumed index of 
refraction of yi =1.33; the imaginary part of ~n was ignored, an 
approximation that should be acceptable in view of the high value 
of €Oo . Values of ot 0 < / would tend to be physically 
uninteresting, because /i would be quite large and the corresponding 
greenhouse effect quite small. On the other hand, values of 
*f > Y would not show much diversity since /3 would uniformly 
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remain about unity, and the single scattering patterns, both in 
the visible and the infrared, would change only slowly. 


TABLE 1 

Single Scattering Parameters 


63o 

a: 




0.99 

0.2 

16 

1 

5 

0.99 

0.4 

32 

2 

2 

0.99 

0.5 

64 

4 

1 


The relevant phase functions in the infrared, demonstrating 
the increase of the forward scattering lobe with increasing oi** , 
are shown in Figure 1, Computations of for, the 

visible, however, cannot be disposed of so easily. It is possible 
to make the necessary computations for the values of oC£ 
listed in Table 1, but the computer time required is rather 
excessive. However, the major difficulty stems from other 
considerations. The extreme nature of the forward scattering 
lobe of ■jpCc&z* <3>) for such large values of cx? requires a 
quite large value of N in ,a least-squares Legendre polynomial 
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fit of p (cos 0) in accordance with equation (3), and the computer 
program developed to solve the corresponding equation of transfer 
becomes unstable for values of N much in excess of N = 17. Re- 
finements in the computational techniques become a matter of 
considerable delicacy beyond this point, and high values of <jJ 0 
tend to compound the difficulties. 

Rather than attack the problem directly it was decided to 
attempt an alternative approach. First, phase functions (71- 1.33) 
were computed from the Mie theory for o<^ = 1, 2, and 4, and then 
fitted by equation (3) for N = 13. Next, the maximum value of 

was sought for which equation (3) would satisfactorily fit 
the rigorously computed phase function. In this way it was found 
that o£ = 6 would yield a phase function that could be- approxi- 
mated by (3) everywhere to within one percent of the true value, 
and much better over the forward scattering lobe; larger values 
of o were found to yield unsatisfactory results. Finally, the 
scattering pattern for very large spheres ( >> l) for a 
refractive index of 1.33 was adopted from Table 21 on page 232 
of van de Hulst's (1957) book. This phase function was calculated 
by van de Hulst from geometrical optics, neglecting diffraction 
effects. 
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The five phase functions mentioned above have been reproduced 
in Figure 2. The solid curves were computed from the rigorous Mie 
theory, except for van de Hulst’s large sphere ( ctf/ =: 00 ) model. 
The dashed curve approximating the last phase function was computed 
from (3) for N = 13; in all other cases the approximate and exact 
phase functions are indistinguishable in the figure. We shall 
see in the next section how these five phase functions are used 
to approximate realistic limb functions and temperature profiles 
for the atmospheric models of interest. 

V. TEMPERATURE PROFILES AND 
THE LAW OF DARKENING 

We shall be interested in solving equations (67) -(68) for 
the temperature profile and (63) for the law of darkening (T-O) 
for each of the models listed in Table 1. (in order to obtain 
absolute temperature determinations in conformity with the solar 
flux at the distance of Venus, we have adopted a solar constant 
ttF 0 = 0.26 watts/cm2jt In addition we will want to obtain solutions 
for each of these models for varying zenith angles of incident 
solar radiation; for this purpose the values /<© = 1.0, 0.4, and 
0.16 have been chosen to be representative of the range of interest. 
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It has already been remarked that the actual computations must 
be restricted to ot**6 excepted). In order to obtain 

approximately correct results for the models listed in Table 1, 
we have introduced the following artifice. 

Consider first the model ( ®t = j 5% = .S'j • fi-l ' ^ 4 s/,0 }; 

for this model c£ = 64 . It is required now to compute the 
temperature profile from (67) and (68) , using instead of = 64 

the values o/* =1, 2, 4, and 6. In the upper, part of Figure 3 
the temperatures thus computed from equations (67) -(68) for 
infrared optical depths of T = 0, .1 ,i;iO, and oo are plotted as a 
function of 1/cX % . It is now assumed that the radiation diff- 
racted around spheres very large compared with the wavelength 
(i.e. - oo ) is contained completely in the forward direction 

( ® = O ) . Physically this is the same thing as saying that the 
diffracted radiation does not contribute to the diffuse radiation 
field. According to the Mie theory, which includes the effects of 
diffraction rigorously, the limiting value of as cx£ — > oo is 
Xe — 2 X$ , where is the geometrical cross-section. One 

half of this limiting value of X& is due to diffraction effects. 
It would appear to be a fairly good approximation in the case of 
large spheres to assume that none of the diffracted radiation is 
lost from the radiation field, since the corresponding cross-section 
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from geometrical optics is just . Hence, a reasonable approxi- 
mation to impose on the model we are concerned with in the 
limiting case of o(£-co is to (1) retain the scattering pattern 
calculated by van de Hulst, (2) reduce by a factor of two, 

and, in order to account for no losses of diffracted radiation, 

(3) reduce ( l — Z3% ) by a factor of two. The resulting sjingle 
scattering parameters become = .S' ‘ o£*r 9 * /3 - »£* }. 

The temperatures resulting from a solution for this model 
for the optical depths of interest are included along the ordinate 

’i 

corresponding to l/o<o * 0 in the upper part of Figure 3. With 
the addition of these last points the curves connecting points of 
equal T can be drawn quite smoothly. In addition, the general 
character of each curve tends to imply that valid temperature 
profiles for any value of o£ / can be successfully extracted 
from the figure,. Of course, the only temperature profile thus 
inferred having any physical significance is the one ( ©</ = 64) 

we are interested in, and is indicated by the arrows in the upper 
part of the figure. 

The lower part of Figure 3 contains a similar family of 
interpolation curves for the model = j £;/= ,2 J e 'fisS }/*<>-,/£} 
For this model the limiting case was calculated using 

values of /3 = 2.5 and = .98. The arrows in the lower part of 
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of the figure refer.this time to the physically meaningful value 
c< 9 = 16. Notice in particular the peculiar nature of the curve 

along T-0 relative to other curves for different values of V . 

The temperature profiles in Figures 4 and 5 refer for the 
most part to the upper and lower families of curves in Figure 3 
respectively. The five left-most solid curves in each of the 
figures are the temperature profiles computed for the five size 
parameters = 1, 2, 4, 6, and co , while the dashed curves are 
the temperature profiles interpolated for each of the two 
physically meaningful values of o<o . It is interesting to note 
that a temperature minimum exists for the cases in Figure 5 some- 
where in the range (.12 T £ ,19). For large zenith angles 

of the sun it appears that, under the proper conditions, the extreme 
upper part of a particulate atmosphere is heated primarily by direct 
solar radiation, while further down the heating is caused primarily 
by the diffuse radiation field. It would seem that the temperature 
minimum occurs where the effects of the two mechanisms are of 
about equal magnitude. The high surface temperatures obtained 
for the sun at the local zenith (Figure 4) are about the same 
as those reported by Clark and Kuz'min (1965) for Venus, and hence 
are quite suggestive that particulate matter may be of some 
significance in the generation and maintenance of a greenhouse 
effect in the atmosphere of Venus. 
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In order to test whether an increase in the forward scatttering 
of visible radiation (resulting in a deeper penetration into the 
atmosphere) could produce a significant increase in the temperatures 
at large optical depths, the case of complete forward scattering 
was considered. Physically this would require that all the 
scattered radiation was scattered into the direction (S> •= O • This 
may be effected mathematically by replacing &C with zero and 
reducing , in accordance with the original value of ££ , by 

the factor (1 - ) . Further considerations reveal that the 

constants /_„<(<#.= 1, ...» n) in the relevant equations [cf. eqs. 

(63) and (67) ] should be set equal to zero, and the' constant Y 0 
identified with unity. The right-most solid curves in Figures 4 
and 5 illustrate the resulting temperature profiles for the 
respective cases in a rather dramatic fashion. Physically these 
cases refer to large spheres having a refractive index of unity in 
the visible (with a small absorption component) . It would appear, 
therefore, that one method of enhancing the greenhouse effect is 
merely to require a small index of refraction in the visible. Of 
course the extreme case of complete forward scattering will result 
in a planetary albedo of zero; hence, as the index of refraction 
in the visible is steadily reduced, the corresponding value*, 
of must be steadily increased in order to maintain a constant 

t 

planetary albedo. 
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Through the method of interpolation described above, tem- 
perature profiles for each of the three models listed in Table 1 
Were computed for solar direction cosines of = 1.0, 0.4, and 
0.16, and the results depicted in Figure 6. It is readily 
apparent that the temperature profile is relatively insensitive 
to the particle size down to an optical depth of *C ** 6 for the 
models adopted, but becomes quite sensitive for T >10. Notice 
that as /8 decreases the effective heating at large T increases, 
as one would expect intuitively. Far more important than /Q is 
the effect of A • The increase of temperature at all optical 
depths with increasing A is due to the corresponding increase 
of incident .solar flux and decrease of planetary albedo. , This 
latter effect is enhanced by the forward scattering nature of the 
individual particles. 

Table 2 below lists the planetary albedos A and effective 
temperatures T„ computed for each of the nine models presented in 

Cl 

Figure 6 in accordance with the relations 

A I F(6) 

A ~ 1 a£ (72) 

and 

£-/ » <r —v 

Ft 0 ) - y T B } (73) 
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where F(0) follows from the flux integral (56). The values of 
I )(i = 1, 7) in (56) are identically zero by virtue 

of the boundary conditions imposed at the top of the atmosphere, 
while the required values of 1(0, /^)(i = 1, 7) were obtained 

by interpolation procedures entirely analogous to those used for 
obtaining the temperature profiles previously discussed [cf. 

Fig. 3]; instead of interpolating along curves of equal T , inter- 
polations were made along curves of equal . These latter 
interpolations weije equally as satisfactory in appearance as the 
former. The results are illustrated in Figure 7 in units of 

temperature. 

# * 


TABLE 2 

Planetary Albedos and Effective Temperatures 
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The albedos in Table 2 appear to be somewhat on the low side 
compared with the observed values for Venus (cf. Kozyrev, 1954), 
while effective temperatures listed in the table average out 
appreciably higher than those observed for Venus (cf. Murray, 
Wildey, and Westphal, 1903) , and the variation of T £ with is 
nothing like that observed. Again, the laws of darkening illus — 
trated in Figure 7 average out to be somewhat more extreme than 
those observed by Murray, et al., while the variation of the laws 
with /a c do not agree with the observations at all. In view of 
the low rotation rate (Carpenter and Goldstein, 1963) and high 
effective nighttime temperatures observed for Venus it is highly 
likely that strong convection and advection are the dominant 
modes of heat transfer in the atmosphere of the planet. Such modes 
of heat transfer would be expected to decrease horizontal tem- 
perature gradients initiated through radiative processes quite 
markedly. However, it is possible that surface temperatures of 
600°K — 7 00 near the equator, and i-S0 0 K ad the. poUs 
exist (Clark and Kuz’min, 1965), implying that large horizontal 
temperature gradients near the surface may well be present. It 
appears very likely in any event that heat transfer mechanisms 
other than radiative ' play a large role in maintaining the actual 
temperature profiles existing in the atmosphere of Venu§. Never- 



theless, the probability of radiative processes being responsible 
in part for generating the gross features of the actual vertical 
temperature distribution of the atmosphere cannot be ignored. 

VI. CONCLUDING REMARKS 

It has been demonstrated that a thick cloud of particulate 
matter suspended in an atmosphere of optically inactive gases is 
capable of supporting a rather large greenhouse effect. In some 
regards this greenhouse mechanism is particularly promising in tfco 
case of Venus. One point that has not been emphasized is con- 
cerned with the small vertical distances required to attain large 
atmospheric optical thicknesses. As an example, let r = 2 x 10” 4 cm 
X e = 2 7 r Z- 2, , and the particle density N Q = 100 particles/cm* 3 . 

A simple calculation in accordance with the rough relation [cf. 
eq. (31) ] 

T-K xf 2. (74) 

‘I 

shows that the height of the atmosphere H© corresponding to an 
optical thickness of V = 100 is about 2© 4 0 km. Hence there is 
no need to ascribe the high surface pressures to the atmosphere of 
Venus that is required by a greenhouse model based solely on the 
infrared blanketing properties of carbon dioxide and water vapor 
(Jastrow and Rasool, 1963). 1 
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It should be noted that Sobolev f s (1963) value of S3* = *989 
was inferred from a much more isotropic phase function in the 
visible than would be expected to be relevant for the atmosphere 
of Venus. This is borne out in the somewhat low planetary albedos 
(computed for S3/ = .99) listed in Table 2 of the present study. 

It is reasonable, therefore, to suspect that the value of for 
Venus is considerably higher. ’.In addition, it is unreasonable to 
neglect gaseous absorption and emission in the infrared in view 
of the substantial amounts of carbon dioxide known to exist in the 
atmosphere of Venus. An inclusion of these effects can only 
serve to reduce the values of used in the present study, since 
absorption in the visible should be relatively minor by contrast. 

A reduction in , in turn, will serve to raise the surface 
temperature without affecting the effective temperature Tg(cf. eq. 
(64)]. On the other hand, increasing will serve to reduce T E> 
although it is not so obvious how the surface temperature will be 
affected. In any event it would appear likely that certain more 

realistic combinations of oS* 0 and /3 than those considered in 

1 f , 

this study would maintain or increase the high surface temperature 
and at the same time reduce Tg, i.e. increase the planetary albedo. 

A quantitative parametric study, more complete than the one 
presented here, as well as. an extension of the general theory to 
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nongrey atmospheres, would be of considerable interest. In 
addition, any effort made to be definitive in describing the 
thermal structure of the atmosphere of Venus without including 
convection as an important mode of heat transfer is probably 
useless. We intend to return to these and related topics at a 
later time. 

The unfailing willingness of Mr. J. D. Barksdale to perform 
the many tasks involved in writing the major computer program and 
providing the required computations is hereby gratefully 
acknowledged. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


- Phase functions for single scattering in the infrared, 

computed for square size parameter distributions over 
the range ( 0 i oi 1 * £ ) . The curves are averages 

of X9 phase functions computed for various complex 
indices of refraction covering the range ( /»/ — 0,01 Ji 

£ ft < A 6 - O.S’O A ) . 

- Phase functions for single scattering in the visible, 

computed for square size parameter distributions over 
the range ( 0 £ ot/ ) . The splid curves are 

the exact Mie calculations for a real index of refrac- 
tion of ^ = 1.33, except for the case of o(% — go 
which was computed from geometrical optics (neglecting 
diffraction). The dashed curve is a 14-term Legendre 
polynomial approximation for the last case. 

- Curves of temperature vs. optical depth used to 
interpolate for two T-T relations for realistic 
size parameter maxima ot% . The arrows for the upper 
family of curves locate the relevant size parameter 

y 

= 64, while the arrows for the lower set locate 
the value = 16. 

- Temperature profiles for the model { = j 

= , £~ j ^ — ‘f j ft— / j ~ 1,0 }. The 

temperature profile corresponding to(,the physically 
significant size parameter maximum = 64 as inter- 
polated from Figure 3 is shown by the dashed curve. 

The right-most solid curve is the temperature profile 
computed for the limiting case of complete forward 
scattering from very large spheres. 
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FIGURE CAPTIONS (Continued) 


Figure 5 


Figure 6 


Figure 7 


A 


- Temperature profiles for the model { CD* zz • 

£%= .2 ; / ' /3 = jr ;//,= ./i }. The 

temperature profile corresponding to the physically 
significant size parameter maximum as inter- 

polated from Figure 3 is shown by the dashed curve. 

The right-most solid curve is the temperature profile 
computed for the limiting case of complete forward 
scattering from very large spheres. Notice the 
transition from a linear to a logarithmic scale at 

r = 1. 

- Temperature profiles for the models listed in Table 1 
as interpolated in accordance with the discussion 
given in the text (cf. also Figure 3). 

- The laws of darkening for the models listed in Table 1 
as interpolated in accordance with the discussion given 
in the text. The corresponding temperature profiles 
are shown in Figure 6. 


39 


I 







©° 







180 190 200 210 220 230 240 250 260 270 280 290 300 310 320 330 340 350 360 370 380 

T (°K) 




250 











